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Center vortices in MCG - finite-size and gauge-copy effects 
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In SU(2) lattice gauge theory in maximal center gauge, we investigate the dependence of center-projected Creutz ratios and 
the vortex density on lattice size and the number of gauge copies. The dependence on the number of copies is rather strong 
on small lattices, but almost disappears on lattices sufficiently large compared to the expected average vortex thickness. The 
center-projected string tension, evaluated on sufficiently large lattices, is in good agreement with the full asymptotic string 
tension, and the vortex density scales according to the second-order asymptotic-freedom formula. 



1. INTRODUCTION 

Center projection in the direct version of maximal 
center gauge (DMCG) [[ljl is a tool for detecting cen- 
ter vortices in SU(A r ) gauge theory. A number of nu- 
merical studies indicate that center vortices are ubiq- 
uitous in the QCD vacuum and give rise to the linear 
confining potential. Recently, however, the validity of 
the procedure based on center projection in MCG has 
been put in doubt by Bornyakov, Komarov, Polikar- 
pov, and Veselov (BKPV) because of the Gribov copy 
problem [||]. Thus we carefully investigate the depen- 
dence of DMCG results on lattice sizes, gauge copies 
and convergence criteria in SU(2) gauge theory. 

In DMCG, the procedure is to maximize 

by an iterative over-relaxation procedure, where 
Tr^t/] is the trace of U in the adjoint representation. 
Let R n denote the value of R after n over-relaxation 
sweeps. When R n is judged to have converged ac- 
cording to a criterion of the form 



then the link variables U^x) are projected onto 
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the nearest center elements Z using Z^{x) = 
signTr^x)]. 

When this gauge-fixing procedure is applied to dif- 
ferent gauge copies of a given lattice configuration, 
it converges not to the unique global maximum of R, 
but to different, local maxima corresponding to differ- 
ent (Gribov) copies. One way to minimize the gauge 
copy dependence is to carry out the over-relaxation 
procedure on a number N copy of random gauge copies 
and perform center projection on the copy with the 
largest value of R. 

Center-projected Wilson loops, Creutz ratios, etc. 
are observables computed from the center-projected 
link variables. It was shown in a number of studies |Jl|], 
that 

• thin vortex excitations of the projected lattice CP- 
vortices") are located roughly in the middle of 
thick center vortices on the unprojected lattice; 

• projected Creutz ratios Xcp(I, I) are close to the 
asymptotic string tension on the unprojected lat- 
tice ("center dominance"); 

• the density of P-vortices, for j3 > 2.3, scales ac- 
cording to asymptotic freedom; 

• removing center vortices (located via the pro- 
jected lattice) from unprojected lattices also re- 
moves confinement and chiral symmetry break- 
ing, and brings the topological charge to zero [^J. 

The question - why should this procedure work at 
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Figure 1. X cp° PV (4, 4) vs. N copy for (3 = 2.5. 

all? - was addressed in ref. There it was shown 
that in the absence of Gribov copies (i.e. if the gauge 
can be fixed to a global maximum of R) center pro- 
jection in maximal center gauge will always locate a 
thin vortex inserted anywhere on the lattice. This was 
dubbed the "vortex-finding property" of MCG. How- 
ever, the vortices in the QCD vacuum are not thin, but 
of finite thickness in physical units, and maximal cen- 
ter gauge is plagued with Gribov copies. 

Thus we will study in this article the sensitivity of 
center-projected Creutz ratios and the P-vortex den- 
sity with respect to: (i) the number N copy of used 
random gauge copies; (ii) the lattice size; and (iii) 
the convergence parameter S in eq. (||). We will here 
only discuss direct maximal center gauge, although 
there exist alternatives like Laplacian center gauge [|J 
(which is free of the Gribov copy problem). 

2. N copy AND LATTICE SIZE DEPENDENCE 

In the original simulations of ref. [|J, only three 
gauge copies were used. Projected Creutz ratios 
Xcp(I,I) were found to be close to the asymptotic 
string tensions reported by Bali et al. [||, for all 
/ > 2. In [§], however, BKPV calculate Creutz ra- 
tios in the range N copy G [1, 20], and extrapolate to 
Ncopy — ► oo by fitting their data to the functional 

form Xc P ooPV {IJ) = XcpilJ) + c(/,/)/ v ^V^. 
They report that projected string tensions, at (3 = 
2.4, 2.5, underestimate the full string tension by about 
20% at N copy = 20, and by 30% in the extrapolation 
to N, 
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Figure 2. Xc P {I, I) for various N CO py, (3 and large L. 

2.3, 2.4, and 16 4 at (3 = 2.5, whereas the data re- 
ported in ref. [|l]j was obtained with lattice sizes 16 4 
at j3 = 2.3, 2.4, and with 22 4 at /3 = 2.5. Therefore 
the BKPV results may be seriously contaminated by 
finite-size effects. To find out, we have repeated the 
center-projection calculation at (3 = 2.3 and (3 = 2.5 
on a variety of lattice sizes, for N CO py G [1, 20]. For 



However, BKPV used lattice sizes 12 4 at (3 = 



the convergence parameter in eq. (|2|), we have used 
S = 2 x 10~ 7 . 

In Fig. [j] we show results for the Creutz ratio 
Xcp c ° Py (4, 4) vs. Ncopy at (3 — 2.5, for lattice sizes 
ranging from 8 4 to 28 4 . As noted by BKPV, there is 
a slow downward trend in Xcp as N CO py increases, but 
this effect is much more pronounced on smaller lat- 
tices than on larger lattices. And although Creutz ra- 
tios on the smaller lattices grossly underestimate the 
full string tension a reported by Bali et al. [Q], the 
data seems to approach a as the lattice size increases. 
These trends in the data are not unique to Xcp(4, 4) at 
f3 = 2.5, but are typical of all of our results. 

In Fig. U we display the projected Creutz ratios 
Xc P copv (I,T) for Ncopy = 5, 10, 15, 20, and for the 
extrapolation N CO p V — > oo on the largest lattices we 
have used: 20 4 at/3 = 2.3, 2.4, and 28 4 at/3 = 2.5. As 
usual in MCG, all the Xcp{l , I) for / > 2 are close to 
a, and these results are not far from our earlier results 
reported in ref. [|lj] . In Fig. |3| we show the average Xcp 
of the projected {N CO p y — > oo) Creutz ratios Xc P {I, I) 



3 



=2.5 



10000 



0.05 
0.04 
0.03 
0.02 
0.01 


-0.01 



Balietal. 




p 











10 



15 20 
L 



30 



1000 



100 













P' 


p 


P 



Figure 3. Xc P (I, I), I = 2-5, for N copy -> oo. 

in the range 7 = 2 — 5. Note the approach of Xcp to a 
as lattice size increases. 

In addition we checked that the numerical results 
are stable with respect to the constant 6 in eq. (fy. 
We calculated \ cv with convergence criteria S = 
1(T 2 , 1CT 3 , 1CT 4 , 2 x l(T 7 at/3 = 2.5 (24 4 lattice). 
For the weakest criterion S = 10~ 2 , Creutz ratios 
come out too high, but the results for the two smallest 
values of S are fairly consistent, indicating that these 
numbers are not far from the 8 — > limit. 

3. VORTEX DENSITY AND THICKNESS 

The lattice P- vortex density p is the number of P- 
vortex plaquettes (i.e. plaquettes on the projected lat- 
tice with ZZZZ = — 1), divided by the total number 
of plaquettes on the lattice, p is proportional to the av- 
erage area taken up by P-vortices per unit lattice vol- 
ume, and is determined from the center-projected pla- 
quette expectation value viap = 1/2(1 — W cp [l, 1]). 
If this quantity scales as predicted by asymptotic free- 
dom, the expression 



p = p 



102/121 



exp 



(3) 



should be constant in the large (3 limit. As can be seen 
from Fig. there appears to be good evidence for this 
kind of scaling for (3 > 2.2. The vortex densities at 
(3 = 2.3, 2.4, 2.5 are taken from the largest used lat- 
tices and extrapolated to the N copy — » oo limit. The 
other values (with N copy = 3) are taken from our pre- 
vious work. It is interesting that the scaling of p, in the 



Figure 4. Rescaled vortex density p = p/F((3). 

range j3 — 2.3 — 2.5, is even better than the scaling of 
the full asymptotic string tension a in this range. 

In order to avoid large finite size effects, rela- 
tively large lattices are required for center projec- 
tion. We think that the relevant length scale is asso- 
ciated with the vortex thickness. Center vortices are 
surface-like objects in D = 4 dimensions which have 
a finite thickness in physical units. The thick vortex 
surface (or "core") bounds a Dirac 3-volume, which 
represents the region of discontinuity of a singular 
gauge transformation associated with the vortex. In 
ref. [Q] we have explained the vortex-finding property 
of MCG in terms of the global properties of this sin- 
gular gauge transformation. On a small lattice, with 
an extension comparable to the vortex thickness, these 
global aspects of the vortex field may be almost ab- 
sent, and the minimal Dirac 3-volume could be quite 
small. In that case, the argument of ref. [Q] breaks 
down. Thus we expect center projection to be less ef- 
fective at finding vortices on small lattices, leading to 
underestimates of both \cp and p. 

If there is some truth in this explanation, center pro- 
jection works well only for lattices which are large 
compared to the vortex thickness. This thickness can 
be deduced from 

1 . the ratio of "vortex-limited" Wilson loops; 

2. the vortex free energy as a function of lattice size; 

3. the adjoint string -breaking length. 

Vortex-limited Wilson loops are defined in the fol- 
lowing way [|l[|: W n {C) is a Wilson loop evaluated 
on a sub-ensemble of unprojected configurations, se- 
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lected so that precisely n P-vortices, in the corre- 
sponding center-projected configurations, pierce the 
minimal area of the loop. For W\ (C) we can further 
require that the negative P-vortex plaquette lies at (or 
touches) the center point of the rectangular loop. It is 
then expected that W\{C) /W (C) — > -1 in the limit 
where the vortex core is entirely contained within the 
loop. At /3 = 2.3 the vortex appears to almost fit in- 
side a 5 x 5 loop [[[]], which leads to a rough estimate 
of the vortex radius of about 3 lattice spacings. A di- 
ameter of 6 lattice spacings at /3 = 2.3 corresponds to 
a vortex thickness of about one fermi. 

A second estimate is obtained from the recent nu- 
merical calculation of vortex free energy vs. lattice 
size [0] . The vortex free energy is close to zero when 
the lattice extension is greater than the vortex thick- 
ness, and this again gives an estimate for the vortex 
thickness of a little over one fermi. Finally, if confine- 
ment is due to center vortices, then an R x T Wil- 
son loop in the adjoint representation must change 
from a (Casimir scaling) area-law falloff to a (color- 
screening) perimeter-law falloff for R greater than the 
vortex thickness [[|. The adjoint string -breaking dis- 
tance has been measured to be 1.25 fm [g], which is 
roughly consistent with the other two estimates. 

At /3 = 2.5, one fermi corresponds to 12 lattice 
spacings. BKPV used at j3 = 2.5 only a 16 4 lattice 
which may simply be not large enough compared with 
the vortex thickness in order to reliably identify vor- 
tices using center projection. 

The average distance between vortices is given by 
the vortex density. The P-vortex density, discussed 
above, is in fact an overestimate of the actual cen- 
ter vortex density because P-vortices fluctuate within 
the thick vortex core [|l0|] . A more accurate estimate 
of the center vortex density is arrived at by either 
"smoothing" the P-vortex surfaces, or else directly 
from the string tension using / = |(1 — e~ a ). The 
two estimates agree fairly well [p0[], and for /3 = 2.3 
we find / = 0.063. This implies an average distance 
of J -1 / 2 ~ 4 lattice spacings between the centers 
of vortex cores piercing a plane. Since we have al- 
ready estimated the vortex thickness at (3 = 2.3 to 
be about 6 lattice spacings, there must be some over- 
lap between vortex cores. There is nothing in princi- 
ple wrong with that; vortex cores are not impenetra- 
ble objects, and their long-range effects are associated 
with Dirac 3-volumes, rather than the detailed struc- 



ture of the core. These findings indicate that the QCD 
vacuum is more like a liquid of vortices than a dilute 
gas. 

4. CONCLUSIONS 

We have found that center-projected lattices are 
more sensitive to finite size effects than unprojected 
lattices. Precision results for center-projected Creutz 
ratios \cp require lattice sizes which are large com- 
pared to the physical vortex thickness of ~ lfm. The 
Ncopy -dependence for \cp reported by BKPV [[|] is 
greatly reduced as lattice size increases. On the largest 
lattices we have used, and with the extrapolation to 
the N copy — > oo limit, our results lie quite close to 
the asymptotic string tension, and are stable with re- 
spect to the gauge-fixing convergence criterion. For 
the vortex density, we find good evidence for asymp- 
totic scaling^ 

REFERENCES 

1. L. Del Debbio, M. Faber, J. Giedt, J. Greensite, 
S. Olejnfk, Phys. Rev. D58 (1998), 094501. 

2. V. G. Bornyakov, D. A. Komarov, M. I. Polikar- 
pov, A. I. Veselov, JETP Lett. 71 (2000), 231. 

3. P. de Forcrand, M. D'Elia, Phys. Rev. Lett. 82 

(1999) , 4582. 

4. M. Faber, J. Greensite, S. Olejnik, D. Yamada, 
JHEP 12 (1999), 012. 

5. C. Alexandrou, M. D'Elia, P. de Forcrand, Nucl. 
Phys. Proc. Suppl. 83-84 (2000), 437. 

6. G. S. Bali, C. Schlichter, K. Schilling, Phys. Rev. 
D51 (1995), 5165. 

7. T. G. Kovacs, E. T. Tomboulis, Phys. Rev. Lett. 
85 (2000), 704. 

8. M. Faber, J. Greensite, S. Olejnik, Phys. Rev. D57 
(1998), 2603. 

9. P. de Forcrand, O. Philipsen, Phys. Lett. B475 

(2000) , 280. 

10. R. Bertie, M. Faber, J. Greensite, S. Olejmk, J. 
High Energy Phys. 03 (1999), 019. 

11. V. G. Bornyakov, D. A. Komarov, M. I. Polikar- 
pov, hep-lat/0009035 (2000). 

2 After the conference, a new paper of Bornyakov et al. Jll| ] ap- 
peared. Using the method of simulated annealing, they find even 
higher maxima of R, eq. flip, and the center-projected string tension 
smaller than the physical one. It seems clear that a modification of 
the gauge-fixing procedure and/or condition is required. 



